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This is the first of two papers to study a new emergent symmetry which connects orbitally excited 
heavy baryons to the ground states in the combined heavy quark and large iV c limit. The existence 
of this symmetry is shown in a model-independent way, and different possible realizations of the 
symmetry are discussed. It is also proved that this emergent symmetry commutes with the large 
N c spin-flavor symmetry. 
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Quantum Chromodynamics is now almost universally accepted as the theory which governs strong interaction. This 
theory has been repeatedly tested against experiment, with great success. Due to its non-abelian nature, however, 
the QCD coupling gets strong at low energy, and the dynamics become nonperturbative and intractable. As a result, 
much of our quantitative understanding of low energy hadron properties are based on symmetry considerations. The 
most notable of these schemes is chiral perturbation theory, which is based on the fact that, when the light quark 
masses m q — > 0, the QCD Lagrangian is invariant under the chiral symmetry group SU(n/)iX SU(n/)ij. In the real 
world, the light masses are not zero; nevertheless chiral symmetry survives as an approximate symmetry of QCD, with 
symmetry breaking terms of order p/ A x or m v /A x , where m w is the pion mass, p is the scale of external probes, and 
A x the chiral symmetry breaking scale. Despite being just an approximate symmetry, chiral symmetry nevertheless 
provides strong constraints on low energy pion dynamics. 

Important insights into some states in QCD comes from emergent symmetries which are not symmetries (not even 
approximate symmetries) of the QCD Lagrangian, but emerge as symmetries of the states in the Hilbert space of an 
effective theory obtained by taking certain limits. A famous example of such emergent symmetries is the heavy quark 
symmetry [[f]-[| for heavy hadron states containing a single heavy quark with mass niQ ^> Aqcd- Heavy quark spin 
symmetry ensures that states related by a heavy quark spin flip, like (B, B*) and (£;,, E£), are degenerate. Moreover, 
Oh, heavy quark flavor symmetry implies that the brown mucks (i.e., the light degrees of freedom) of heavy hadrons are 
insensitive to the mass or the flavor of the heavy quark. This guarantees that the B — > £)W and A;, — > A c semileptonic 
form factor (which are usually referred as Isgur-Wise form factors) are normalized to unity at the point of zero recoil, 
1 where the initial and final hadrons have the same velocity. Such absolute normalizations of form factors have profound 
implications in experimental determination of the CKM matrix element V c &. The combined heavy quark spin- flavor 
symmetry is described by the symmetry group SU(2ng) where no is the number of heavy flavors, and this symmetry 
is broken by corrections proportional to powers of Aqcd / mQ ■ 

Note that heavy quark symmetry is not a symmetry of the QCD Lagrangian; if it were, the B and D mesons, 
related by heavy quark flavor symmetry, would be degenerate. However, if our interest is restricted to states with 
a single heavy quark, one can perform a spacetime-dependent phase redefinition such that the heavy quark mass 
and spin drop out of the Lagrangian. In other words, while heavy quark symmetry is not a symmetry of the QCD 
Lagrangian, it is a symmetry of the Lagrangian of heavy quark effective theory, which describes only states with a 
single heavy quark. 

Another well-known emergent symmetry is the light quark spin-flavor symmetry for baryons in the large N c limit. 
The large N c limit was first studied by 't Hooft for mesons and was subsequently extended to baryons by Witten [|| . 
They studied how QCD amplitudes involving various numbers of mesons and baryons scale with the number of color 
N c when N c is large. Gervais and Sakita [§,^0) realized that, for large N c baryons, the spin symmetry SU(2) and flavor 
symmetry SU(n/) (where nf is the number of light flavors) are combined and enlarged into the spin-flavor symmetry 
group SU(2ny). (This spin-flavor symmetry was systematically re-analyzed by other groups; see Refs. |Il]^15|].) It was 
shown that the low-lying baryon spectrum in the large N c limit contains a tower of states with (J, J) — (0, 0), (1, 1), . . . 
when N c is even, and (i, i), (|, §),••• when N c is odd. In the latter case one can identify the (5, i) state as the 
nucleon, and the (|, §) as the A(1232) resonance. Moreover, it can be shown that the splittings between these tower 
states are of order 1/N C . As a result, when N c — ► 00, the nucleon, A and all the other states in the tower collapse 
into degeneracy, signifying the emergence of the symmetry SU(2riy). Similarly, in the heavy baryon (baryon with 

a single heavy quark) sector, this spin-flavor symmetry decrees that the Eg -Aq splitting vanishes in the large N c 
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limit. Again, note that this light quark spin- flavor symmetry is an emergent symmetry in the sense that it is not a 
symmetry of the QCD Lagrangian, but only a symmetry of the QCD Hamiltonian of states with unit baryon number. 

We have recently reported H] a new emergent symmetry of QCD which emerges in the heavy baryon (baryon with 
a single heavy quark) sector in the combined heavy quark and large N c limit. This contracted U(4) symmetry (or 
more generally U(d+ 1) in a theory with d spatial dimensions) connects the ground state heavy baryon to some of its 
orbitally excited states, which become degenerate with the ground state as tuq — > oo and N c — > oo. As a result, static 
properties such as the axial current couplings and the moments of the weak form factors of these orbitally excited 
states can be related to their counterparts of the ground state. While many of these results have been discussed 
before in the literature in the context of particular models, they were first presented as model- independent symmetry 
predictions in Ref. |l4f| . 

After the publication of Rcf. |i4[ , we realized that this contracted U(4) can be further enlarged into a contracted 
0(8) symmetry (or more generally 0(2d + 2) in a theory with d spatial dimensions). Moreover, this "symmetric 
realization" is only one of two possible realizations of the emergent symmetry. In the "symmetry broken realization" , 
the symmetry is broken down to contracted 0(4) x 0(4). This paper is the first of two papers where we will report 
these and other new progresses on these emergent symmetries in the combined heavy quark and large N c limits, as 
well as discuss the results reported in Ref. ]l4|] in more details. This paper will focus on the formalism and the different 
realizations of the emergent symmetry, while phenomenological applications and corrections to the symmetry will be 
discussed in the next paper 

This paper is organized as follow: In Sec. II, we will briefly review the bound state picture, a class of models which 
exhibits this same contracted 0(2d + 2) symmetry and which motivates our studies. While the bound state picture 
is not logically related to QCD, it provides a simple physical picture of the origins of this new symmetry. The bound 
state picture treats the heavy baryon as a bound state of an ordinary baryon and a heavy meson and thus is a model. 
However, emergent symmetries are often first recognized in models. For example, heavy quark symmetry was first 
discovered in quark models, while the large N c spin-flavor symmetry was first realized in the Skyrme model. Hence it 
is useful to first consider a model which embodies the correct symmetry to get a feeling of the physical picture before 
launching a formal discussion of the symmetry in QCD language. 

After examining the logical foundation of the bound state picture in Sec. Ill, we begin the main task of this 
paper and study the emergent symmetry in the context of QCD. In Sec. IV, we discuss the relative sizes of different 
contributions to the QCD Hamiltonian. Then in Sec. V we introduce the kinematic variables of the bound state 
picture in the context of QCD, and show that many conclusions of the bound state picture can be justified in the 
model independent manner. The generators of the emergent symmetry will be formally introduced in Sec. VI, and in 
Sec. VII and VIII we will focus on the "symmetric realization" and show that in this case the QCD Hamiltonian is 
that of a three-dimensional simple harmonic oscillator by considering multiple commutation relations in the combined 
heavy quark and large N c limit. Following this is a short discussion in Sec. IX, while Sec. X and XI will discuss the 
"symmetry broken realization" and the inclusion of spin and isospin effects. Then the paper concludes with a short 
preview of the companion paper , which is under preparation and will discuss phenomenological issues and higher 
order corrections to the symmetry predictions. 



II. THE BOUND STATE PICTURE OF A HEAVY BARYON 



The bound state picture fll6|-p4|] regards a heavy baryon as the bound state of a heavy meson and a light baryon (a 
baryon without any valence heavy quarks); the latter often treated as a chiral soliton, i.e., a topologically nontrivial 
configuration of the classical meson fields. In particular, the lightest charmed baryon A c is regarded as the bound 
state of the heavy mesons D or D* (which are degenerate in the heavy quark limit) and a nucleon. In the following, 
we will focus on the model described in Refs. |l6|-|i~8|] , which will be referred to as the simple bound state model as 
it is the simplest model with correct behaviors in the heavy quark and large N c limit. However, we emphasize that 
one can make a similar analysis on other versions of the bound state picture, and the symmetry properties should be 
qualitatively the same as long as these models are consistent with heavy quark symmetry and obey the usual large 
N c counting rules. 

In QCD, a heavy baryon is a complicated bound state, with the quarks interacting through strongly coupled gauge 
dynamics, and with quark-antiquark pairs popping in and out of the vacuum, etc. — a highly intractable problem. 
The bound state picture replaces it (in an ad hoc manner) with the much-simpler problem of a two-body bound state. 
Moreover, the problem further simplifies in the heavy quark limit, where the heavy meson becomes infinitely massive, 
and the large N c limit, where the nucleon mass mjv grows like N c . For concreteness, we will adopt the prescription 
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(only for this section) that the heavy quark limit is taken before the large -ZV C limit. [] Taking the heavy quark limit 
first, the reduced mass of the two-body system [i ~ mjv ~ N c — > oo in this combined heavy quark-large N c limit. As 
a result, the kinetic term, which is suppressed by vanishes, and the wave function does not spread but is instead 
localized at the bottom of the potential. (When /j, — ► oo, the absolute square of the wave function will be a Dirac 
delta distribution at the bottom of the potential.) Consequently, a small attraction between the heavy meson and the 
nucleon is sufficient to ensure the existence of a bound state. 

What is the potential V(x) between a heavy meson and a nucleon? Without resorting to models, we do not know 
much about the potential except the fact that, by the usual large N c counting rules V(x) is of order N®. However, 
let us assume that V(x) has a global minimum at the origin, i.e., the heavy meson sits on the top of (the center 
of) the nucleon. In this case, the wave function will be highly localized around the origin, and the potential can be 
approximated by V(x) = Vq + ^kx 2 , which includes only the first two terms in the Taylor expansion of V(x). For the 
origin to be a global minimum, we need Vq < and k > 0. In this case, when the bound state is the ground state 
of the simple harmonic oscillator, it is a Aq. On the other hand, excited states in the simple harmonic oscillator are 
orbitally excited heavy baryons. With explicit wave functions, coupling constants and form factors for transitions 
between different states can be calculated in a straightforward manner. 

However, it remains to be seen whether the assumptions that Vq < and k > are justified. In the simple 
bound state model [^6|-|l8|, the nucleon is described as a chiral soliton, i.e., a topologically nontrivial classical pion 
configuration. By using a chiral Lagrangian which is truncated to the leading order, the potential energy of a heavy 
meson in the presence of a classical background pion field can be calculated. It turns out that indeed Vq < and 
K > 0, and the assumptions are verified in this particular model. As a result, one can identify different heavy baryons 
with eigenstates in a simple harmonic potential. However, this result is clearly model dependent. 

As mentioned before, the bound state picture possesses an emergent symmetry which relates the ground state to 
the excited states. This can be seen by making the crucial observation that the excitation energy ui — \J k/ fx is small, 
where fi is the reduced mass of the bound state and k is the spring constant of the simple harmonic potential. By 
first taking the heavy quark limit, /i = m^v (mass of the light baryon) scales like N c . On the other hand, since the 
binding potential V(x) is of order N®, the spring constant k — being its Taylor coefficient — is generically also of 

order N®. Hence ui scales like iV c and vanishes in the large N c limit. This implies that when N c — > oo, the whole 
tower of excited states becomes degenerate with the ground state — a signature of an emergent symmetry. 

What is the symmetry group of this emergent symmetry then? It has to contain, as a subgroup, the symmetry 
group of a three-dimensional simple harmonic oscillator, namely U(3) generated by = a\a,j (i,j — 1, 2, 3) where 
a,j is the annihilation operator in the j'-th direction. These T^ 's satisfy the U(3) commutation relations. 

[Tij,Tki] = SkjTu — SuTkj. (2.1) 

Note that this U(3) group contains the rotational SO(3) subgroup, generated by Li = —itijkTjk with [Li, Lj] — ie^^Lk- 
When N c — ► oo and the excited states become degenerate with the ground state, the annihilation and creation 
operators aj and a\ (i,j = 1, 2, 3) also become generators of the emergent symmetry. The additional commutation 
relations are 

[a>j,Tki] — Skjai, [a\,T M ] = ~5ua\, [aj,a\] = 6^1, (2.2) 

where 1 is the identity operator. These sixteen generators {Tij, ai, at., 1} form the minimal spectrum generating algebra 
of a three-dimensional harmonic oscillator, i.e., the smallest algebra which contains the symmetry group U(3) and 
connects all eigenstates of a three-dimensional simple harmonic osc illat or. It is related to the usual U(4) algebra, 
generated by (i,j = 1, 2, 3, 4) satisfying commutation relations ( |2.l| ) by the following limiting procedure: 

a, = lira T 4j /R, aj = lim T l4 /R, 1 = lim T 44 /R 2 . (2.3) 

R — >oo R — >oo R. — >oc 

Such a limiting procedure is called a group contraction, and hence the group generated by {Tij, ai,a\, 1} is called a 
contracted U(4) group. ^ 



1 In the real world, the heavy meson masses tnj ~ 5 GeV, mo ~ 1.8 GeV while the nucleon mass mjv ~ 1 GeV. So as far as 
heavy baryon kinematics is concerned, the real world is closer to the heavy quark limit than the large -/V c limit, justifying our 
ordering of the limits. 

2 This contracted U(4) is different from the contracted SU(4) group of the light quark spin-flavor symmetry Jn|. 
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The contracted U(4) minimal spectrum generating algebra can be enlarged to contain the extra operators Sj. 



and ST 
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= 1, 2, 3). The new commutation relations are 

[S^, Ski] — [Sij,ai] = 0, [Sij,a k ] — a^jk + ajSik, [Sij,Tja] — SuSjk + SjAk, 
[S^, S kl ] = TikSji + TjkSa + TuSjk + TjiSik, 



(2.4) 



and the commutation relations involving Sjj can be obtained through Hermitian conjugation. As a result, these 
28 generators {Sij, Sj-, Ty, ai, aJ k , 1} form a closed operator algebra, which is actually a contracted 0(8) algebra. 
This contracted 0(8) is generated by the creation and annihilation operators, all possible bilinears, as well as the 
identity operator. Note that one cannot further enlarge this operator algebra by including trilinears in ctj and aj; 
the commutator of two trilinears, for example, will be a quadralinear, and the algebra will not close (or will contain 
an infinite number of generators). As a result, this contracted 0(8) algebra is called maximal spectrum generating 
algebra of the three-dimensional simple harmonic oscillator. Again, as N c — > oo and u> — ► 0, the excited states become 
degenerate with the ground state and the contracted 0(8) become the symmetry group of the bound state picture. 
The relationship between all the algebraic structures discussed above is summarized in the following chain: 
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maximal spectrum 
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as u) — » 0. 



(2.5) 



We have shown that the contracted 0(8) is a symmetry of the bound state picture. In the more general case of 
a bound state picture with d spatial dimensions, it is clear that the symmetry is described by a similarly contracted 
0(2g?+ 2) group with a contracted U(d+ 1) subgroup. While we have been focusing on the simple bound state picture, 
this symmetry is actually a generic feature of all variants of the bound state picture as long as the models exhibit 
heavy quark symmetry as toq — > oo, and obey the large N c scaling rules as N c — > oo. However, it is not obvious 
that the physical picture is intuitively reasonable. This will be addressed in the next section. 



III. THE FOUNDATION OF THE BOUND STATE PICTURE 



Questions may be raised about the foundation of the bound state picture on several different levels. On the technical 
level, one may question the description of a nucleon as a classical pion distribution in the simple bound state model. 
Because we have infinitely many species of mesons in the large N c limit jjJU , there is no reason why all other mesons 
besides pions should be ignored. This question can be resolved in part by including more light mesons in the model. 
This is the motivation behind Ref . 0] , where the effects of the p and oj vector mesons are included, leading to results 
which are numerically improved at the expense of more parameters and much more mathematical complexities. Since 
we are interested in the generic features of the bound state picture, we will only remark that including extra meson 
states does not change the physics qualitatively. However, in the large N c limit there is an infinite number of mesons, 
and each meson has infinitely many coupling constants. 

A more serious technical issue of concern is the modeling of the interaction between the heavy meson and the 
classical light meson fields (which make the nucleon). In the simple bound state picture, the heavy mesons interact 



3 A note on the literature: as far as the authors can discern, among all the literature on the bound state picture for heavy 
baryon, only the simple bound state picture Jl6|-[l8[ makes the explicit statement that the binding potential is simple harmonic 
in the combined heavy quark and large N c limit. In none of these works on the bound state picture was the spectrum 
generating algebras discussed, nor was the observation that in the combined limit they become the symmetry group of an 
emergent symmetry. Both of these points were first explicitly raised in Ref. J^]. On the other hand, the appearance of an 
emergent symmetry in the bound state picture does not depend on the details on the model, as long as the model embodies 
the heavy quark and large N c symmetry. Consequently, the emergent symmetry is an implicit feature of all viable bound state 
models. 
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with the classical background pion configuration through a truncated chiral Lagrangian, which is the most general 
interaction Lagrangian which respects chiral symmetry truncated to the leading order of the chiral expansion p/A x , 
where p is the pion momentum. While it is justifiable to use this truncated Lagrangian in low momentum pion 
processes (where p <C A x ), there is little justification for such truncation here, as the chiral soliton contains pionic 
modes over a wide range of momentum, and in general p/A x is not a small expansion parameter. As a result, the 
truncated Lagrangian is an ad hoc description of the interaction between the heavy meson and the classical pion fields. 
The situation is even worse in models where the p and a vector mesons are included. Since the interactions involving 
these vector mesons are not well-constrained by symmetry (unlike pionic interactions, which are severely constrained 
by chiral symmetry), their interactions with heavy mesons are only described by phenomenological Lagrangians of an 
entirely ad hoc nature. 

However, these technological issues are not fundamental and do not alter the conceptual issues about the bound 
state picture. For example, it seems likely that, as far as the emergent symmetry is concerned, the description of the 
light baryon as a chiral soliton is not essential. The essence of the bound state picture is that the heavy baryon can 
be regarded as a bound state with potential V(x) ~ N° and reduced mass \i ~ N c . The details of the interaction are 
inessential as far as the symmetry is concerned. This naturally leads us to ask the question whether one can recast 
the analysis in such a form that chiral solitons are not invoked. As we will see below, the answer to this question is 
affirmative. 

A more severe conceptual criticism of the simple bound state model is the use of point particle quantum mechanics 
when both the heavy meson and the nucleon are extended objects. Assuming that the bound state picture is reasonable, 

the mean square relative displacement of the nucleon from the heavy meson is 3/(2/xw) ~ N c , which vanishes as 
N c — > oo, while the size of the nucleon has a smooth non-zero large N c limit. Hence the heavy meson will be jiggling 
well inside the nucleon near its center, and it is not obvious that point particle quantum mechanics is applicable. 

Lastly, the connection of the bound state picture to QCD is obscure. To address this philosophical concern, one 
can only try to reproduce the emergent symmetry directly from QCD. This is the purpose of both our previous paper 
and this paper. We will see that in a model-independent way, one can show that this contracted 0(8) symmetry 
is not only a symmetry of the bound state picture, but in fact a symmetry of QCD. 



IV. DISSECTING THE QCD HAMILTONIAN FOR HEAVY BARYONS 

Due to the conservation of baryon number and heavy quark number (in the heavy quark limit), it is legitimate to 
restrict our attention to the heavy baryon Hilbert space, i.e., the subspace with both heavy quark number and baryon 
number equal to unity. In the combined heavy quark and large N c limit, this subspace is well-defined. We introduce 
the small power counting parameter, A, to quantify the deviation from the combined limit. It is defined as: 

A-AOER 1 ( 4.i) 

with the ratio NcAQcn/m-Q arbitrary. In other words, both 1/toq and 1/N C corrections are of order A, while order 
A 2 corrections include those scale like 1/toq, 1/wqJV c and etc. 

In the heavy baryon Hilbert space, it is useful to decompose the QCD Hamiltonian H in the following way: 

H = Hq + He, where Hq — niQ + Hq and He = mjv + He. (4.2) 

The heavy quark part of the QCD Hamiltonian Hq contains the heavy quark mass uiq , as well as the heavy quark 
kinetic and interaction terms denoted by 'H.q. Since Hq involves only a single quark (namely the heavy quark), it at 
most scales like N® ~ A . (In contrast, both ttlq and tun are large and of order A" 1 .) As we are only interested in 
states with a single heavy quark, by performing a Foldy-Wouthuysen transformation Hq can be expanded in powers 
of 1/ ttlq in heavy quark effective theory: 

H Q = gA° + g^- + 0(m Q % (4.3) 

2m<g 2m Q w 

where Pq is the three-dimensional heavy quark momentum: 

Pq= I d 3 xh(x)Dh(x), (4.4) 
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with h(x) being the heavy quark field in heavy quark effective theory, and D is the three-dimensional covariant 
derivative. Note that the chromomagnetic term Sq ■ B is suppressed by the heavy quark mass rriQ. While the 
definitions of toq and Pq are ambiguous since the heavy quark mass is not uniquely defined, these ambiguities are of 
order unity (A ) while the heavy mass and mo ment um are typically large (toq, Pq ~ A -1 ). Q As a result, the relative 



ambiguities are small and one can rewrite Eq. (4.3) as 



Hq = gA° + A. - + 0(A 2 ). (4.5) 

2uiq 2niQ 

Similarly, the light part He contains the nucleon mass (which is proportional to N c ) and He, which represents 
the change in the energy of the brown muck (i.e., the light degrees of freedom of the heavy baryon) when one of the 
light quarks is replaced by a heavy quark. We cannot write down a simple expression for He as we have done for Hq , 
but it is easy to see that it scales like N®. The reasoning is as follows: the interaction energy between any two quarks 
is of order N~ x by the standard large iV c counting rules Since the replacing of a light quark with a heavy quark 

in a baryon breaks N c — 1 light quark-light quark interactions and replaces them with N c — 1 light quark heavy quark 
interactions, we have in the large N c limit, 

He ~ (number of interactions modified) x (change of energy in each interaction) ~ N c X N' 1 ~ N° ~ A . (4.6) 

The light Hamiltonian He contains all the dynamics of the brown muck as well as its interaction with the heavy 
quark. In general, it can depend not only on its position x and momentum p relative to the heavy quark, but also all 
kinds of internal degrees of freedom which correspond to different modes of excitation. In comparison, the Hamiltonian 
H hs of a two-particle bound state, in general, can be decomposed in the following form: 

H hS — Kkin + Hpot + Hcxc, (4.7) 

where 7ikin is a kinetic term which depends only on p, H po t is a potential term which only depends on x, and H e ^ c 
represents possible intern al ex citations and commutes with both x and p. The issue becomes whether He can be 
recast in this form of Eq. fl4.7| ) . This is the question which we will be attempting to answer in the next four sections. 

V. KINEMATICS AND THE KINETIC ENERGY 

In the previous section, we have decomposed the QCD Hamiltonian, H, into a heavy part Hq and a light part 
He- Our next step will be to perform similar decompositions for the kinematic variables; namely, the momentum 
and position operators. We will reproduce the two-body kinematics of the bound state picture using QCD operators. 
Recall that our aim is to demonstrate the existence of an emergent symmetry in QCD itself. In order to achieve this 
goal in a model-independent manner, one cannot simply assume the kinematic variables in the bound state picture 
are well defined. Instead we need to construct these kinematic variables from QCD operators without reference to 
any model. 

While the total momentum of any given heavy baryon system P is a well-defined quantity, in general there is 
no unambiguous way to separate the momentum into a heavy quark contribution and a brown muck contribution. 



However, in the heavy quark limit, the heavy quark momentum Pq in Eq. (4.4) is a well-defined QCD operator (up 
to corrections of relative order A), and one can define the brown muck momentum Pe as P — Pq. Lastly, the QCD 
based position operators of the whole system X, of the heavy quark Xq, and of the brown muck Xe are defined as 
the conjugate operators of the respective momentum operators: 

[Xj,P k ] = -iSij, [X Qj ,P Qk ] = -idij, [Xtj,Pt k ] = -iSij, [xj,Pk] = -iSij, (5.1) 

where Xj is the j-th component of X, etc.. In the last equality, x is defined as the relative position operator Xe — Xq, 
and p is its conjugate operator. The relationship between these eight operators (four momenta and four positions) 
are summarized in the following diagram: 



4 Both the heavy quark mass ttiq and the heavy quark velocity v are well defined up to order A ambiguities. For a discussion 
of these ambiguities and their theoretical implications, see Refs. Eq,BQ|. 
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P = P e + P Q p 

t I I I. (5-2) 

X Xi — Xq = x 

where vertical arrows represent conjugations. Of course these look just like the analogous relations in the bound state 
picture, but recall that the point of introducing these operators is to see whether the bound state picture dynamics 
can be reproduced directly from operators in QCD with no model-dependent assumptions. 

By construction, the heavy quark momentum and position operators commute with the brown muck counterparts. 

[X Qj ,X tk ] = [X Q .,P lk ] = [P Qj ,X tk ] = [P Qj ,P lk ] = 0. (5.3) 

The center-of-mass position X is an unknown linear combination of Xq and Xg. Similarly, the relative momentum p 
is an unknown linear combination of Pq and Pg. The operators are defined in such a way that the relative kinematic 
variables commute with the center-of-mass counterparts. 

[X, , x k ] = [X 3 , Pk ] = [P 3 , x k ] = [Pj , p k ) = 0, (5.4) 

which in turn implies the following linear relations: 

X = aX e + $X Q , p = $P e -aP Q , (5.5) 

with a and [3 being operators which commute with all the momentum and position operators and satisfy a + (3 = 1, 
the identity operator. We emphasize that all of these operators arc defined from the QCD operators P and Pq through 
linear combinations and conjugations. As a result, all of them are QCD operators. 

What are the operators a and 0? To answer this question, one needs to look at the dynamics of the system and 
study the QCD Hamiltonian TL. In particular, we will study the commutators of TL with the position operators. 
First, consider the commutator [Xj,TL\. 

[X j ,H]=iX j = i5L, (5.6) 

where the second equality is from Poincare invariance. ^ Note that the Pj/TL is well defined as TL commutes with 
Pj. Moreover, since TL = M + TL, where A 
relative correction of order A. As a result, 

[Xj,H] = iPj/M + C(A 2 ), (5.7) 

and consequently we have the following double commutators: 

[X k , [X^H]} = -5 jk /M, [x k , [X jt H]] = 0, (5.8) 

where all 0(A~ 2 ) corrections are dropped. Note that the first equality is exactly what one expects if one starts with 
the nonrelativistic Hamiltonian, P 2 /2M. As a result, we say that the kinetic mass of the heavy baryon, defined to 
be the reciprocal of the double commutator with the position operator, is M = rriQ + in the combined heavy 
quark and large N c limit, with possible corrections of order unity. These double commutators will be important in 
the determinations of a and (3. 

Next, let us study the commutator [Xq.^TL]. Out of the four contributions to the QCD Hamiltonian TL in Eq. ( |4.2|) , 

rriQ and mjy are c- numbers, and the light operator Hi commutes with Xq.. Thus, 

[X Q H] = [X Q Hq] = [X Q gA" + 9%-^] = + 0(A 2 ), (5.9) 

J 3 3 2mn 2mn ran 



Pj. Moreover, since 7i = M + H, where M — itlq + tojv ~ A 1 while Ti ~ A , one can replace 1/Tt with 1/M with 



5 For the reader who does not find the above equality obvious, recall that TC commutes with Pj (this is one of the defining 
commutation relations of the Poincare group), and hence the rest mass M, defined by TL — ^2 P 3 + M , is Poincare invariant. 
Moreover, being the energy of the whole system at Pj — 0, M is given by the QCD Hamiltonian TL in Eq. (|l.2|); as a result M — 
m Q + m N up to corrections of order A . As a result, [Xj,TL] = i-jp-TL = i-jp-(Y, k P l + m2 ) 1/2 = iP J '(J2k P k + M 2 )- 1/2 = i^L. 
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as A and B are light operators and commute with Xq. On the other hand, by ignoring the 0(\ 2 ) corrections, one 
has the following double commutators: 

[X Qk , [X Qj ,H]] = -Sj k /mQ, [X lk , [X Qp H]] = [X Qk , [X £j ,H]] = 0, (5.10) 

where the Jacobi identity has been used to show the vanishing of the last double commutator. Note that the first 
equality states that, in the heavy quark limit, the kinetic mass of the heavy quark is tuq. 

Si nce X and x are linear combinations of Xq and Xg (cf. Eq. ( |5.5D ), one can recast the double commutators in 
Eq. (5.8) in the following form: 

- 5 jk /M = [X k , [X J -,W]]= a 2 [X lkl [X tj ,H]] + 2 [X Qk , [X Qj ,H]}, 

0= [x k ,[X j} 7{\] = &[X lk ,[X tj ,7{\]- 0[X Qk ,[X Q .,7{\], (5.11) 

where the cross terms vanish by the second formula in Eq. ( 5.1 0| ) . We know [XQ k , [Xq^,TC\] from Eq. (5.10) but do 
not know [Xi k , [Xij, TL\] at this stage. Canceling the latter, we end up with 

- 5 jk /M = (a + $)$[X Qk , [X Qp H]} = -$6 3k /m Q , (5.12) 

where the relation a + (3 = 1 has been used. Finally, this implies 

f3 = m Q /M, a = m N /M. (5.13) 

Thus the operators a and (3 turn out to be c- numbers in the combined heavy quark and large N c limit. 

What do these values of a and (3 tell us about the QCD Hamiltonian, Til First, it is now straightforward to show 
that the kinetic mass of the brown muck is tun- 

[X lk , [X tj ,H]\ = -S jk /m N + 0{\ 2 ). (5.14) 



This, together with Eq. (5.1C), implies the following decomposition of 7i: 

H = m Q + m N + H Q xm + Hi,Un + Hpot, (5.15) 
where in the combined heavy quark and large N c limit, 

n QM n = ^- + 0(X 2 ), 7^, kin = T ^- + 0(X 2 ), (5.16) 
Zing ZrriN 

where the second equality comes from Eq. ( |5.14 ). The potential energy term 7Y po t ~ A does not depend o n an y 
momentum operators (but can and does depend on the position operators; see the following section). Using Eqs. (5.13), 
these two kinetic terms can be recast in terms of the center-of-mass and relative momenta. 

P 2 f 

"Hkin = T^Q.kin + W^kin = + — , (5-17) 

where M = mq + is the total mass and fi = mq-m^ / (mg + m^) will be referred as the reduced mass of the 
system. (Note that both M and are of order A -1 .) The reduced mass /i can be interpreted as the kinetic mass of 
the relative coordinate. Indeed, from the obtained value of a and /3, one can verify that 

[x k ,[ Xj ,H}} = -6 3k /fi + 0(X 2 ) (5.18) 



is in agreement with Eq. ( 5.17 ). In other words, in the combined heavy quark and large N c limit, the kinetic terms 
of the heavy baryon system are those of a nonrelativistic bound state of two point particles with uiq and m n . 

One may wonder what is the point of this whole exercise of reproducing elementary two-particle quantum mechanics, 
with apparently no new results. However, remember there is no a priori justification of treating a heavy baryon as 
the bound state of two point particles. In particular, the brown muck is not a point particle; it has a substantial 
size and complicated internal structures with the possibility of excitations. A formalism such as the bound state 
picture which treats the brown muck as though it were a point particle requires justification from QCD. Our goal is 
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to demonstrate the existence of an emergent symmetry in QCD itself (not merely in a model), we have to work with 
operators P and Pq, which are QCD operators, and construct out of them all other kinematic operators. That our a 
and (3 are identical with that in a nonrelativistic point particle treatment means that we have succeeded in providing 
a justification of the latter treatment in the combined hea vy quark and large N c limit. 

The apparently trivial double commutator in Eq. ( |5.14 ), which states that the kinetic term of the brown muck is 
nonrelativistic in the combined limit, is in fact not completely trivial. In the presence of a heavy quark, the kinetic 
mass of a composite object like the brown muck may depend on the relative position of the brown muck to the heavy 
quark. What we find instead is a constant kinetic mass in the combined limit. Physically this reflects the fact 
that the system in question is weakly bound; the interaction term TL, which is of order A , is much smaller than the 
masses of the constituents which are of order A -1 . In a weakly bound state, the kinetic mass of the whole system 
is the sum of the kinetic masses of the constituents. Since the kinetic mass of the whole system M and that of the 
heavy quark itiq are position independent, the kinetic mass of the brown muck is also independent of its position. 



VI. THE GENERATORS OF THE EMERGENT SYMMETRY 

In the previous section, we have analyzed the kinetic terms of the QCD Hamiltonian, TL, by studying its double 
commutators with the position operators. One may also want to analyze the potential term 7i po t by studying the 
double commutators of TL with the momentum operators. Unfortunately, this strategy actually provides very limited 
information. 

As Poincare invariance demands that [Pj,TL] = 0, it immediately follows that [Pgj,TL] = — [Pq^TL]. Moreover, 

\ Pj ,'H\ = 0[Pe j ,'H\-&[P Qj ,T{\ = (& + 0)[P ij ,7{\ = [P tj ,H], (6.1) 

where again we have used a + /3 = 1. This reflects the simple observation that TL po t can depend on the relative 
position x but not the center-of-mass position of the whole system X. 

Recall that we have made the following decomposition: TL — ttiq +rriN + 7Ykm + ^pot, with the last two terms both 
being of order unity or less. Both ttiq and mjv are c- numbers and commute with any operator, and it is easy to see 



that 7ikin commutes with all momentum operators from its expression in Eq. (5.17), which does not depend on any 
of the position operators. As a result, 7i po t is the only term which may have non- vanishing commutators with the 
momentum operators. The commutators of interest are [pj, TL po t}, \pk, [Pj^pot]], \jPk, \jPj, [Pi, Wpot]]], etc. We can say 
very little about these multiple commutators except that they are all (at most) of the same order as 7i po t, i-e., of 
order unity. 

Of particular interest is the double commutator, whose significance lies in the following definition of the spring 
operator k: 

KS jk = -\p k , \pj,H]] = -\p k , [Pj.Hpot]]- (6-2) 

Let |G) be the ground state of the QCD Hamiltonian TL in the heavy baryon Hilbert space, and Eq be its mass, 
satisfying (TL — Eq)\G) = 0. The spring constant k is then defined as (G\k\G). Since [] 

k = -fa, [ Pj ,TL]} = -[p 3 , [ Pj , {TL - E )]} = 2 Pj (TL - E ) Pj - PjPj {TL - E ) - (TL - E ) PjPj , (6.3) 
it is easy to see that n is positive by inserting a complete set of states {!".)}• 

k = 2(G\ Pj (H - E ) Pj \G) = 2J2 \(n\Pj\G)\ 2 (E n - E ) > 0. (6.4) 

n 

Now we arc ready to define the generators of the emergent symmetry; namely, the creation and annihilation 
operators. 



In the following two equations, the repeated index j is not being summed over. 
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where \i = mQrriN / {rriQ + m^v) is the reduced mass and ui — y/n/fx will be referred to as the natural frequency of the 
heavy baryon. With k of order unity and // ~ A -1 in the combined heavy quark and large iV c limit, us vanishes — an 
important result which does not depend on the order in which the two limits are taken. If the heavy quark limit is 

— 1/2 

taken first, /i ~ N c and ui ~ N c ' — > as N c — * oo. On the other hand, if the large N c limit is taken first, /i ~ rriQ 

and u> ~ jrq 1 ' 2 — > as tuq — ► oo. 

The natural frequency u> plays a central role in heavy baryon dynamics in the combined heavy quark and large N c 
limit. As we will see below, ui is the coefficient of the only term in the QCD Hamiltonian which breaks the emergent 
symmetry at leading order of A. As a result, all the physical properties (masses, coupling constants, form factors, 
etc.) of the low-lying baryons can be expressed in terms of ui. Alternatively, if one can determine the value of to by 
measuring some physical observable which depends onw {e.g., the mass of the first excited heavy baryon), then one 
can predict the values of many other physical observables. 



VII. CONSTRAINTS ON THE QCD HAMILTONIAN 

Let us recall that our analysis is based on A counting of quantities describing heavy baryon dynamics in the combined 



limit. The reduced mass can be expressed as the double commutator in Eq. ( 5. IS ) 



5jk/fi = ~[xk, [xj,H]\ ~ A, (7.1a) 

where higher order terms in A are dropped. Similarly, the spring constant k is the ground state expectation value of 
the double commutator in Eq. (6.2): 

5 jk k = -\p k ,\ Pj ,H}]~ \°, K = (G\k\G). (7.1b) 

As a result, the natural frequency uj = y^JJi ~ A 1 / 2 . This is a notable feature: recall that the expansion in A 
embodies the expansions in both Aqcd/^q and 1/N C . We rarely encounter situations in which fractional powers of 
Aqcd/wq or 1/N C arise for direct physical observables. Here, however, we have found that powers like A 1 / 2 do arise 
naturally. In fact, we will see that the natural expansion parameter will be A 1 / 2 instead of A. This ultimately reflects 
the interplay of two heavy scales (namely ttlq and tun) 



The double commutation relations in Eqs. (7.1) constrain the possible forms of the QCD Hamiltonian TL. Note, 
however, that both double commutation relations are satisfied by replacing TL with 7isHOi the Hamiltonian of the 
bound state picture, which is just the simple harmonic oscillator. 

Hsho = | + — - T =«"»}• a,. (7.2) 

Moreover, the contracted 0(8) symmetry mentioned above is precisely the maximal spectrum generating algebra 
of TisHO and becomes an emergent symmetry as ui — > as A — > in the combined limit. On the other hand, to 
demonstrate that this contracted 0(8) is a symmetry of QCD in this limit, one needs to show that the generators of 
the contracted 0(8) commute with the QCD Hamiltonian H, or equivalently, show that 

H = WsHO + ^exc + • ■ ■ , (7.3) 

where 7i oxc commutes with a and a* in the combined limit, i.e., [aj,7i cxc ] = [aj,7i oxc ] = 0, and represents the 
possibility of internal excitations of the brown muck. On the other hand, the ellipses are possible corrections to the 
simple harmonic Hamiltonian TLsao and should be negligible in comparison to TLsko m the combined limit. In other 
words, we need to show that, relative to 7isHO; these correction terms are suppressed by powers of A, and hence can 
be dropped in the combined limit. 

Before we embark the power counting of TL, we need to clarify the counting of powers for the kinematic variables. 
For example, consider the simple harmonic Hamiltonian TLsno = wo' • a, which apparently is of order A 1 / 2 as u> ~ A 1 / 2 . 

However, TLsho can also be written as ^- + ^| where the kinetic term, being suppressed by l//i, is apparently 

of order A, while the potential term, with coefficient k ~ A , is apparently of order unity. The origins of this apparent 
discrepancy lie in the relationship between the operators (x,p) and (a,a>). 



^-(a + at), j?=-i^( S -rf). (7.4) 
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Since fju ~ A -1 / 2 , one cannot simultaneously set x, p, a and eft to the same order in A. 

In the following, we will make the prescription that a, ~ A , which implies x ~ A 1 / 4 and p ~ A -1 / 4 , and show 
that it is self-consistent. As we will see later in this paper, this prescription will lead to the "symmetric realization" of 
the emergent symmetry with a contracted 0(8) symmetry group. We can make the following a posteriori justification 
for this "symmetric prescription" . Instead of studying the power counting of the operators x, p, a and eft , one can 
study the power counting of the matrix elements of x, p, a and eft between the low-lying states of TL. Power counting 
on matrix elements, which are c-numbers, is free of the aforementioned ambiguities. If the low-lying states of TL are 
simple harmonic states as suggested by the bound state picture, then the matrix elements of a and a) are indeed of 
order unity, while the matrix elements of x and p are not. As the subsequent discussion verifies this picture, we will 
have justified our prescription a posteriori. 

Now we are ready to show that in the combined limit TL has the form presented in Eq. ( |7.3| ) . We will achieve this 
in three steps. In the remainder of this section, we will verify the fact that all possible triple commutators of TL with 
a and aft vanish in the combined limit. More specifically, we will show that these triple commutators go to zero more 
quickly than HshOj which scales like A 1 / 2 . Then in the following section, we will show how the vanishings of these 
commutators imply that TL can be at most bilinear in a and a) : 

TL = Caft -a + Da-a + D^aft -a) +H CXC , (7.5) 

where TL™ r commutes with both a and a T . Lastly, C and D can be determined from the double commutation relations 
in Eqs. (|Tl|). 

The relevant triple commutators of TL with a and eft are the following operators: 

t<°> = K [ aj , [a k ,TL]}}, tW= [aj, [a h [a k ,TL]}], 

t^ = [4, [at, [a k ,TL]]], t&= [aj, [a], [a£,H]]]. (7.6) 

Using the Jacobi identity and the commutators of aj and at, it is easy to show that the values of these t^ triple 
commutators do not depend on the ordering of the Oj's and aj-'s. For example, 

[al [a h [a k ,0]]] = -K-, [[a k ,0],a\}] - [[a k ,0], [aj,^]] = [a,-, [aj, [a k ,0]]] 

= - [a,, [a k , [0,a\]}} - [ aj , [O, [aj,a fe ]]] = [a h [a k , [a\,0}}}. (7.7) 

Since t^ = -(i (3) ) f and *W = -(i (2) ) t , it suffices to show that t^ = t^ = in the combined limit. 

Since a and are linear combinations of x and p, and can be expressed as linear combinations of the triple 
commutators of TL with x and p: 

t (0) = 2 -3/2 (T (3) + 3zT (2) _ 3T (1) _ iT (0) ); t (l) = 2 -3/2 (T (3) + fT (2) + T (l) + ^(0)^ (? g) 

where 

T<°> = (^r 3 ^ [ Pj , [ Pk ,TL]}}, T«= ( M w)- 1/2 [Pi, fe, [a;fc,W]]], 

= (j^f/ifa, [ Xj , [x k ,TL]]], T&= (^) 3 / 2 [^, fo, [**,«]]]. (7.9) 

Again note that the values of these triple commutators do not depend on the ordering of the Xj's and p/s. We 
have shown in Sec. VI that the triple p commutator is at most of order unity, and hence ~ C(A 3 / 4 ). All of the 
other three triple commutators are also small as [x k ,TL] — ip k //J. + 0(A 2 ). The first term gets killed by the following 
commutations and does not contribute to the triple commutators. So ~ A 9 / 4 , ~ A 7 / 4 and ~ A 5 / 4 — all 
vanish faster than T^°\ As a result, both t^ and t^ vanish at least as fast as ~ A 3 / 4 , and are negligible when 
compared to TLsko ~ A 1 / 2 in the combined limit. 

VIII. THE QCD HAMILTONIAN IN THE COMBINED LIMIT 

In the previous section, we have verified that all triple commutators t^ vanish faster than TLsw ~ A 1 / 2 in the 
combined limit. In other words, if we only keep terms up to those of order A 1 / 2 , all these triple commutators are zero. 

[a l7 [a jt [a k ,TL]}] = [aj, [a,, [a k ,TL]}} = [aj, [at, [a k ,TL]}} = [aj, [at, [a+,W]]] < 0(\-^ 4 ). (8.1) 
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Now the vanishings of these triple commutators at order A 1 / 2 imply that to this order the double commutators 
commute with a and eft: 

[a j ,[alH}] = [a],[a k ,H]] = -C5 jk , [oj, [a£, H}] = 2D6 jk , [ aj ,[a ki H]]=2ft8 Jh , (8.2) 

where the operators C and D commute with both a and aft. Note that C is hermitian while D is in general non- 
hermitian. The above relations can be recast as: 

[a,j , B] = 0, [a] , B] = 0, where B = [a\.,H] + Ca\ + 2Da k . (8.3) 

However, parity invariance of TL implies that there does not exist any parity odd operator which commutes with both 
a and eft . So B vanishes and 

[a\,H] + Ca{ + 2Da k = 0. (8.4a) 

Similarly, one has 

[a k ,H]-Ca k -2&a\ = 0. (8.4b) 



Notice that, while in Eqs. (3.2) we expressed the double commutators of TL in term of C and D, in Eqs. (8.4) we 
managed to express the single commutators of TL in term of C and D. We can make one more step and express TL 
itself in terms of C and D by noting that Eqs. (f^J) imply 



i]=0, [at,l]=0, where A = H - Ccft ■ a + Da ■ a + D ] ct ] ■ a\ (8.5) 



where A commutes with both a and a> . Lastly, by renaming A as TL CXC , we recover Eq. (7.5): 



H = Ca t ■a + bd-a + D ] a ] ■ a f +H CXC . (8.6) 

We have succeeded in showing that the vanishings of the triple commutations of 7i imply that TL is at most a 
bilinear in a and eft. Q While the above derivation looks rather complicated, the essence of the statement is very 
intuitive: if TL contains, for instance, trilinear terms in a and aft , the triple commutators will read off the coefficients 
of these trilinear terms and hence will not vanish. What we have achieved, through the derivation above, is to realize 
this intuition in a rigorous manner. 

We have shown that, up to order A 1 / 2 , TL is a bilinear in a and eft. Now we will make the final step in this 
demonstration and determine the form of the operators C and D. Clearly if C = lo and D — 0, TL will simply be the 
sum of the simple harmonic Hamiltonian Hsho and a possible excitation term 7i cxc which commutes with a and eft . 
Our goal is to show that these are indeed the values of C and D. 

One can re-express TL in terms of x and p: 

W= — p H x - — — (x ■ p + p- x) +TL C xc, (8.7) 

where t) + = D + £) f and ib- = D — £)t. Then C, D+ and L>_ can be deduced by using the double commutation 



relations Eqs. (7.1) 



\ Pj ,[x k ,H}]= D- =0, 

[xj, [x k ,TL]] = - Sjk/fj, => C-D + =uj, 
\pj, \p k ,TC]] = — §j k k => C + D + = k/fnj = Lok/n. 



The last two equalities lead to: 



7 This is actually a special case of the following general result. Let O be an arbitrary operator with vanishing m-th multiple 
commutators with a and eft . Then one can prove by induction that O is at most (m — l)-linear in a and eft . 
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C = oj(k/k + l)/2, D + =uj(k/k-1)/2, 



(8.8) 



where the spring constant k is the ground state expectation value of the spring operator k introduced in Eq. ( p.2j ). 
As a result, when acting on the heavy baryon ground state, k/k = 1, which in turn implies C — to and D + = 0. 
For a general heavy baryon state (not necessarily the ground state), however, k is not identical to its ground state 
expectation value n. However, it is true not merely for the ground state, but also for states in the ground state band, 
which is the subspace spanned by states of the form (aj,)™ x (o-t) n ' \G). We therefore conclude that, in the ground 

state band, C = uj, D + = 0, and TL has the simple harmonic form: 

n =^ + ^2 2~ + Wcxc + ° (A) = uS S + Wcxc + ° (A) - (8J) 

This Hamiltonian clearly reduces to a three-dimensional simple harmonic oscillator with reduced mass li, spring 
constant n, and hence natural frequency u. The maximal spectrum generating algebra of this Hamiltonian is a 
contracted 0(8) algebra, and in the combined limit when u ~ A 1 / 2 — * 0, the contracted 0(8) becomes an emergent 
symmetry. Again, we emphasize that this emergent symmetry is a symmetry of QCD in the heavy baryon sector — 
not only that of the bound state picture or any other models. 



IX. REVIEW AND DISCUSSION 



In summary, we have demonstrated that the contracted 0(8) symmetry seen in the bound state picture is in fact a 
symmetry of QCD. Near the combined limit there exists a band of low-lying heavy baryons, labeled by (n x ,n y ,n z ), 
the number of excitation quanta in the x, y and z directions (or alternatively (N, L, L z ), where N — n x + n y + n z 
is the total number of excitation quanta, L the orbital angular momentum and L z its z-component). For each state 
the excitation energy is {n x + n y + n z )uj = Nlu. As A — * in the combined limit, uj — ► and the entire band become 
degenerate. Our discussion can be generalized in a straightforward manner to the case with d spatial dimensions, 
with an emergent contracted 0(2d + 2) symmetry group. 

While the demonstration was rather long, the basic idea is very simple. We started by constructing the kinematic 
variables Xj and pj, which are not a priori well defined in QCD (Sec. V). Since our goal is to study the symmetry in 
QCD itself, we cannot merely assume that the kinematic variables are well defined (as in models), but need to show 
that they are legitimate QCD operators. After defining these kinematic variables in QCD, we showed that they behave 
like the quantum mechanical kinematic operators for bound states of two point particles. Thus they can be used in 
a "quantum mechanics-like" framework to describe the dynamics in the heavy baryon sector. With the creation and 
annihilation operators constructed from Xj and pj (Sec. VI), we show, by considering triple commutators (Sec. VII), 
that the QCD Hamiltonian TL is a bilinear of these creation and annihilation operators up to a certain order in the A 
expansion (Sec. VIII). Lastly, the "coefficients" (which are formally operators) of the bilinear terms in TL are fixed by 
considering double commutators. 

It may seem strange that a symmetry of QCD is only applicable to a certain subspace (namely, the ground state 
band of the heavy baryon subspace) of the whole QCD Hilbert space. This, however, is a typical feature for emergent 
symmetries. The heavy quark symmetry is only applicable to states containing a heavy quark and the light 

quark spin-flavor symmetry in the large N c limit is only relevant for baryonic states [pHL3|. While the symmetries of 
the QCD Lagrangian is applicable to all states in the QCD Hilbert space, emergent symmetries are not symmetries 
of the QCD Lagrangian and may be applicable only to particular subspaces. 

We will end this section by briefly returning to the bound state picture and discuss some intricate issues on its 
relationship to our formalism. We noted in Sec. Ill that there were conceptual problems associated with the bound 
state picture, particularly in regard to treating the brown muck as though it were a point particle. The possible 
problem was that the characteristic size of the brown muck distribution is L 2 ~ Aqcd ~ A while the bound state 
wave function had a typical size of ~ (fiuj) 1 ^ 2 ~ A 1 / 2 which is characteristically narrower in position space. In spite 
of this concept we have shown that as far as kinematics and symmetries are concerned, the point particle description 
correctly reproduces the QCD result. It turns out that the comparison between the characteristic size of the brown 
muck and the scale of the wave function is not the appropriate comparison. In fact, there are three distance scales 
in this problem: If the brown muck is to be approximated by a point particle, the point particle should be located at 
the center-of-mass of the brown muck in order to reproduce the correct kinematics. According to Ref. 0, the brown 
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muck can be studied under the Hartree picture, which becomes exact in the large N c limit. [] In the Hartree picture, 
one can easily see that the center-of-mass of the brown muck is much better localized than each individual quark. The 
size of the wave function of each individual quark L q is comparable to the size of the whole brown muck, which is of 
order unity. On the other hand, the fluctuation of the position of ccntcr-of-mass Lqm of the brown muck is smaller 
by a factor of VN C by the central limit theorem. As a result, L 2 CM ~ L 2 /N c ~ A, which is much smaller than the 
typical spread of the wave function. The different scales form the following hierarchy: 

L 2 ~A° » L^-A 1 / 2 » Lq M ~ A. (9.1) 

In other words, the bound state wave function cannot resolve the fluctuation of the center-of-mass of the brown muck. 
This provides an intuitive explanation why, despite its huge size, the brown muck can be approximated by a point 
particle without drastically altering the kinematics and the symmetries of the system. 



X. THE SYMMETRY BROKEN REALIZATION 



Recall that our demonstration of the simple harmonic form of the QCD Hamiltonian H depends on the symmetric 
prescription for the A power counting introduced in Sec. VII, that the creation and annihilation operators should be 
counted as order unity in TL. This prescription is in turn a posteriori justified by the fact that the matrix element of 
dj and at between states in the ground state band are indeed of order unity. While this confirms the self-consistency 
of this symmetric prescription of A counting rules, it does not preclude the possible existence of other self-consistent 
counting schemes. In this section, we will briefly describe other possible realizations of this emergent symmetry 

To gain physical insight, it is useful to think in terms of the bound state picture. Clearly our simple harmonic 
oscillator obtained by assuming that a, eft ~ A° implies that the ground state expectation value (G\x 2 \G) ~ A 1 / 2 
vanishes in the combined limit. This reflects that, as the reduced mass fi — > oo, the center of the brown muck gets 
more and more localized around the heavy quark. This is the scenario where the origin of the relative position space 
x = minimizes the potential energy globally. However, it does not need to be the case. The potential may have 
a "mexican hat" shape with the global minimum at r = |aT| = r > 0, where by naturalness r ~ 1/Aqcd ~ A . 
In such a case, in the combined limit the relative wave function will be a shell sharply peaked around r = r$. As 
a result, (G\x 2 \G) — r 2 , ~ A and the symmetric prescription is clearly inapplicable. Instead this "mexican hat" 
scenario corresponds to a different realization of the emergent symmetry, which hereinafter will be referred to as the 
"symmetry broken realization" . 

We will sketch how one may study the emergent symmetry in this "mexican hat" scenario, where there are two 
modes of low-energy excitations. Firstly, there are orbital excitations along the bottom of the potential well at r = ro, 
described by the Hamiltonian Hl = L 2 /2T. The moment of inertia X — fir^ ~ A -1 in the combined limit. As a 
result, as X 0, Hl ^ and the whole tower of orbitally excited states collapses into degeneracy. The symmetry 
group of "Hl with finite moment of inertia X is the rotational group 0(3), and the spectrum generating algebra is 
contracted 0(4) (also known as E$, the three dimensional Euclidean group; cf. Sec. III. 4 of Ref. P7J). Secondly, there 
are radial excitation around r — ro, which can be studied through a formalism similar to what we constructed in 
previous sections to study the symmetric realization. It turns out that the radial excitations are also simple harmonic 
in the combined limit, but only as in a one-dimensional oscillator again with u> ~ X 1 / 2 . The spectrum generating 
algebra in this case is again a contracted 0(4) (generated by a, a) , a 2 , aJ , a^a and 1). Near the combined limit 
where A is small, the rotational excitation energies ~ 1/2X ~ A are much smaller than the radial excitation energies 
~ uj ~ A 1 / 2 . Hence the spectrum consists of a tower of equally spaced simple harmonic levels with splitting u>, with 
each level further split into a tower of rotor states with splitting ~ 1/2X. As A — * 0, all these orbitally and radially 
excited states become degenerate with the ground state, and the spectrum generating algebra contracted 0(4) x 0(4) 
becomes the emergent symmetry group in this symmetry broken realization. 

This contracted 0(4) x 0(4) group in the symmetry broken realization is a subgroup of the contracted 0(8) in the 
symmetric realization. Actually this is very reminiscent of spontaneous symmetry breaking in field theory. Note that 
orbital excitations are light in the combined limit. As a result, one can construct a wave function sharply peaked at 



Actually Ref. |g| used the Hartree picture to study a baryon, not the brown muck of a heavy baryon. However, the difference 
between a brown muck of a heavy baryon, with N c — 1 light quarks, and a light baryon with iV c light quarks, becomes negligible 
in the large N c limit. (Of course, the brown muck is different from a baryon as the former is not a color singlet, but the extra 
color charge is neutralized by the heavy quark.) 
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x = y = 0, z = ro which is degenerate with the ground state (and hence is itself also a legitimate ground state) as 
A — > 0. Such a ground state breaks the contracted 0(8) in the symmetric realization to the contracted 0(4) x 0(4) 
in the asymmetric realization. This explains the terminologies, "symmetric" and "symmetry broken" realizations. 

While this symmetry broken realization of the emergent symmetry may not be as aesthetically appealing as the 
symmetric counterpart, we emphasize that it is a viable logical possibility and there is no theoretical justification of a 
priori rejecting this possibility. However, these two realizations are phenomenologically distinguishable, as least when 
A — ► 0. In the symmetric realization the excitation energy of the second excited state is twice that of the first excited 
state, where in the symmetry broken realization the ratio is 3. In the real world, the first excited charmed baryon is 
around 330 MeV heavier than the ground state and about 200 MeV beneath the D-N dissociation threshold. If future 
experiments find the second excited charmed baryon beneath this dissociation threshold, one would be very tempted 
to rule out the asymmetric realization. 

In this section, we have compared the two possible realizations of the emergent symmetry: the symmetric realization 
when the potential is globally minimized at the origin, and the symmetry broken realization when the global minimum 
of the potential is away from the origin. (Actually it is logically possible that there is more than one global minimum — 
one at the origin, while the other is not. This scenario, however, is so extremely unnatural and requires such fine-tuning 
that we will not consider it further in this paper.) While the preceding discussion relied on the bound state picture, 
one can rephrase it in QCD language in a manner analogous to our analysis of the symmetric realization represented 
in previous sections. The symmetry broken realization has the interesting feature of light states of excitation energies 
~ A, which originate from the orbital revolution around the bottom of the "mexican hat" potential. But let us be 
reminded that there are other possible rotational modes for a baryon which has nothing to do with orbital revolution. 
For example, one can rotate the brown muck itself (not around the heavy quark) in space or in isospace, which 
quantum mechanically correspond to spin and isospin excitations. The moment of inertia is of order N c , which leads 
to a whole tower of rotor states which are degenerate in the large N c limit. This is the well-known large N c spin-flavor 
symmetry [p|-[l3"l, relating A(1232) to the nucleon and Eg"* to Aq. These rotational modes of the brown muck have 
little to do with its relative motion relative to the heavy quark. As a result, intuitively we expect this large N c 
spin-flavor symmetry to commute with the emergent symmetry (in either realization). It will be the goal of the next 
section to demonstrate that this intuition is indeed correct. 



XI. INCLUSION OF THE SPIN AND ISOSPIN EFFECTS 

So far we have neglected the spin and isospin of the heavy baryon, which consists of a single heavy quark and N c — 1 
valence light quarks. For concreteness we will only consider the cases where N c is an odd number, so there is an even 
number of valence light quarks in a heavy baryon. In QCD with two light flavors, each light quark is isospin-^, and as 
a result the brown muck can be of isospin 1 = 0, 1, ... (N c — l)/2. The isospin symmetry is described by an SU(2)j 
group, generated by 

. w c -i 

I a = d 3 x Ik^Qk, a = l,2,3, (11.1) 
J k=i 

where the summation is over all the valence light quarks. Similarly, light quarks are spin-i fermions; a brown muck 
with N c — 1 light quarks without any orbital angular momentum between them can be of spin S& = 0, 1, . . . (N c — 1)/2. 
(Note that in the heavy quark limit the heavy quark spin Sq decouples from the rest of the system. As a result, the 
brown muck spin Si is conserved and is a good quantum number.) The brown muck spin symmetry is also described 
by an SU(2)s { group, generated by 

r N c -1 

S}= d 3 x vWlk, i = l,2,3. (11.2) 

^ k=l 

Both isospin and brown muck spin symmetries are symmetries of the QCD Lagrangian (the latter only in the heavy 
quark limit), and their generators I a and S\ satisfy these commutation relations: 



9 To make such a conclusion, however, one has to check if the corrections higher order in A are small. Unfortunately, such 
corrections are likely to be substantial. 
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[p , I b ] = ie abc I c , [S}, Sl\ = te^Sj, [I a , SI] = 0. (11.3) 

It was realized in 1993 by several different collaborations [fil]-[l3| that for large iV c baryons, the separate isospin 
and brown quark spin symmetries, described by the product group SU(2)/X SU(2)g f , get combined and enlarged into 
an emergent spin-flavor symmetry. This spin-flavor symmetry is described by a contracted SU(4) group, generated 
by {X al , I a , S}}, satisfying the following commutation relations: 

[X ai ,I b ] = ie abc X ci , [X ai , Sj] = ie ijk X ak , [X al ,X hj ] =0. (11.4) 
with X m being the axial current couplings: 

. Nc-l 

I a! = ^glrVft, 0,1 = 1,2,3. (11.5) 
fc=i 

Since all generators commute with K 2 , where K is the vector sum of / and Si, all the states with the same K 2 

are degenerate. Of particular interest are the K — states, for which I = Si and hence (/, Si) = (0, 0), (1,1), . . . , 

((N c — 1)/2, (N c — 1)/2). Phenomenologically one can identify the (0,0) state as Aq and the (1,1) state as Eg , and the 

Sq^-Aq splitting is l/N c suppressed. These analyses have been extended to orbitally excited baryons in Refs. [p8|-pl|, 
although orbitally excited baryons containing charm or bottom quarks were only briefly discussed. 

Note that X ai , I a and Sj act only on the light quarks. In other words, they represent the internal degrees of 
freedom of the brown muck, while leaving the heavy quark alone. On the other hand, the creation and annihilation 
operators which generate the contracted 0(8) symmetry represent collective excitations; i.e., all the light quarks are 
excited in a correlated manner relative to the heavy quark. Intuitively, one expects these two modes of excitations 
to be independent of each other in the large N c limit, and hence the light quark spin-flavor symmetry group should 
commute with the contracted 0(8) group, q_] 

To show that this piece of intuition is correct, we will prove the following general result: Let Ji be a local operator 
which acts only on the brown muck, i.e., [Xqj, Ji] = [PQji Ji] — 0. Then the forward matrix element of [ctj, Ji] and 

[aj, Ji] all vanish up to order A. For forward matrix element we mean the kinematic condition that there exists an 
inertial frame such that both the initial and final states are at rest (carrying zero total momentum). 
To see that this is true, note that 



[a j7 Ji] = 



and the vanishing of the x and p commutators on the right-hand side imply the vanishing of the a and commutators 
on the left-hand side. It is straightforward to show from the definitions of the kinematic variables that 

x = \{X-X Q ), p = pP-P Q , (11.7) 



where a and f3 are the operator- valued coefficients introduced in Eq. ( |5.5| ). Note that, since a is a c- number up to 
order A, 1/d is well defined up to the same order. Since [Xq-, Ji] = [PQj, Ji] = 0, one has 

[x 3 ;Jt] = \[X 3 , Ji], [ P] , Ji] = p[Pj,Ji). (11.8) 
a 

We have expressed the commutators of the relative kinematic variables x and p in terms of their center-of-mass 
counterparts X and P. However, since Ji is a local operator without any intrinsic momentum scale, [Xj, Ji] — 0. On 
the other hand, [Pj , Ji] in general does not vanish as the local operator is translated in position space by P. The 
forward matrix element, however, trivially vanishes. As a result, both terms in Eqs. ( 11.6| ) are zero, and the proof is 
completed. 



10 Do not confuse the axial current couplings X a ' with Xj, the center-of-mass position! 



11 We note in passing that models based on the bound state picture have often implicitly assumed that spin-flavor degrees of 
freedom commute with the orbital ones. In the simple bound state model Jl^-p^], for example, this assumption was built in 
by using the same profile function to describe the light baryon as a chiral soliton. 
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Two comments are in place here. First, by saying that the commutators [ctj,J(] and [aj, Jg] vanish, we actually 
mean the more precise statement that these commutators are smaller than the typical matrix element of Jg by at least 
an order in A — the order at which it is no longer justifiable to treat a and /3 as c-numbers. Second, by reversing 
the role of the brown muck and the heavy quark, it is straightforward to show that, for a local operator Jq which 
acts only on the heavy quark, i.e., [Xgj, Jq] — [Pgj, Jq] — 0, the forward matrix element of [aj, Jq] and [aj, Jq] also 
vanish up to order A. Both of these results reflect the intuitive statement that any excitation which involves only 
one of the constituents but does not transfer any momentum will not change the relative motion. These results will 
be useful when we consider the heavy baryon matrix elements of pionic current or weak interaction currents in the 
companion paper E3j. 

Returning to our discussion of spin and isospin effects, since the spin-flavor symmetry generators X al act only on 
the light degrees of freedom, they commute with the creation and annihilation operators implying that all of the 
states below are degenerate in the combined limit. The states are labeled by the quantum numbers (N, L, Se, Je), 
where TV is the number of excitation quanta in the simple harmonic oscillator, L is the orbital angular momentum, 
Se is the spin of the brown muck (which is always equal to the isospin / for the K = states which we are working 
on) . Jg = Sg + L is the total (including both orbital and spin) angular momentum of the brown muck. The total spin 
of the whole heavy baryon J = Sq + Jg is the vectorial sum of Jg with the heavy quark spin Sq . Q 



N = 2 I, = 0,2 



N = 1 L= 1 



N = L = 



ly Q2 



A 
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(11.9) 



x az 



The A sector 
J = St = 



The E sector 
/ = S> = 1 



That X a% and a}- commute means that one can construct, say, the Egj state through (^(AT^Aq) = a^sW, or the 

opposite order X al (a^h.Q) — X^AqI, and both constructions give the same state Eg]. 

Recall that we have shown in Eq. (3.9) that the QCD Hamiltonian in the heavy baryon sector can be written as 
the sum of the simple harmonic part HshOi and the internal excitation Hamiltonian Ti cyic , with possible corrections 
of order A. ^ Among the different contributions to 7i cxc is Tij = crl 2 , the Hamiltonian describing the low-lying 
spin- flavor excitations. The large N c spin-flavor symmetry implies a ~ N~ ~ A. As a result, 7ij ~ A and one can 
move it from 7i GXC to 7Ysho anc ^ rewrite Eq. fl8.9p as follows: 



ft = ^SHO + Kxc + O(A), WsHO = ftsHo +ft/, (11.10) 
and 7ioxc IS f ne Hamiltonian for other internal excitation modes excluding the low- lying spin-flavor excitations. 



Here we have many different angular momenta adding in different ways. We will clarify their meanings by comparing 
a heavy baryon to a multi-electronic atom. If the heavy quark is the analogy of the heavy nucleus and the electron cloud 
corresponds to the brown muck, then Sg corresponds to the electronic spin S, L corresponds to the orbital angular momentum 
C, and Jg — Sg + L translates into J = S + C. The heavy quark spin Sq is the counterpart of the nuclear spin I, and the total 
spin of the heavy baryon given by J = Sq + Jg is the analogy of the ^-spin, T — I + J . 
13 We are working with the symmetric realization here. The case for asymmetric realization can be studied in a similar way. 
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The Hamiltonian W SHO describes the low-lying heavy baryon spectroscopy up to order A corrections. (Note that 
these corrections are of the same order as TLi.) Under Tignc eacn of t ne simple harmonic bound states of Hsho is 
split by TLi into a whole tower of states with different /. The masses of the heavy baryon states under TL'^o are 

to = Ag + Nuj + al 2 + 0(X), (H-H) 

where N = n x + n y + n z is the total number of excitation quanta, and we are adopting the customary abuse of 
notation that the symbol of a state also represent its mass, e.g., Af, = m\ b — 5624 MeV. This mass relation implies 
that the orbital excitation energies are the same in the A sector (I — 0) and the E sector (I — 1), i.e., 

c s - lo a = (E« - E«) - (Ag - A Q ) = 0{\). (11.12) 

Formally it also implies spin-flavor excitation energies of the N = 1 states are the same as their N = counterparts: 

2<n - 2 CTo = (Eg - Ag) - (Eg - A Q ) = 0(A). (11.13) 

Unfortunately this relation is actually devoid of information, as both <j\ and cfq are of order A, which is the order of 
the leading order corrections. 

Note that the qualitative features of the low- lying heavy baryon spectrum is correctly given by 7i SHO , which 
specifies the low-energy spectroscopy up to order A 1 / 2 . The not-yet-specified correction terms at order A are small 
when compared to 7ig HO and hence can only perturb the spectrum but not change it qualitatively. On the other hand, 
to make quantitative predictions about heavy baryon spectroscopy (or other heavy baryon dynamical properties) at 
order A, one has to determine the explicit forms of the order A correction terms — a task which we will undertake in 
our next paper JT5| . 



XII. CONCLUDING REMARKS 



In this paper, we have shown that there is an emergent symmetry in the heavy baryon Hilbert subspace in the 
combined heavy quark and large N c limits. This emergent symmetry can either be realized as a contracted 0(8) 
or a contracted 0(4) x 0(4), and in either case it relates the orbitally excited states with the ground state. Both 
realizations of this emergent symmetry have interesting spectroscopic implications: in the former case the spectrum 
is that of a three-dimensional simple harmonic oscillator, while in the latter case the spectrum is that of a heavy rigid 
rotor with simple harmonic radial excitations. Moreover, we have also shown that this emergent symmetry commutes 
with the light quark spin-flavor symmetry for large N c baryons. 

While the main purpose of this paper is to discuss the results reported in Ref. JIJ] in a more detailed fashion, there 
are several places where the formalism in this paper differs from the original treatment in Ref. . We list the most 
important differences below for comparison: 

• In Ref. Jh|, we have always taken the heavy quark limit before the large N c limit. In this paper, the combined 
limit is taken by keeping NcAqcd/i^q fixed at an arbitrary value. This is a more general treatment, as it includes 
the possibilities of taking the heavy quark limit both before and after the large N c limit, and a whole range of other 
possible limiting procedures. 

• Since the heavy quark limit is taken first in Ref. [ ^4| , one does not need to distinguish the relative momentum 
p from the brown muck momentum Pi. On the other hand, in the paper we are keeping the ratio N c Aqcd /wq 
arbitrary, and the distinction between p and Pe should not be overlooked. Here we have presented the analysis of the 
kinematics in full generality, while the result of Ref. jl4| can be recovered by setting a = and /3 = 1. 

• We have realized that the emergent symmetry group reported in Ref. H], namely the contracted U(4), is only 
a subgroup of the full symmetry group in the symmetric realization — contracted 0(8). Moreover, we have studied 
the symmetry broken realization for the sake of generality. 

We want to emphasize once more that, even though we seem to be dealing with a quantum mechanical potential 
model, the formalism is actually completely field theoretical, and the kinematic variables are QCD operators. Why 
can one reduce a field theoretical problem into a quantum mechanical framework? The answer lies in the separation of 
scales: both constituents, namely the heavy quark and the brown muck, have mass of order A -1 , while the interaction 
is only of order unity and hence is very weak compared to the mass scale represented by the reduced mass /i. The 
separation of scales makes it possible to make an expansion in powers of A and have an effective field theory which 
includes only the lowest excitation modes, which in this case is the motion of the brown muck relative to the heavy 
quark (and the possible spin- flavor excitations). This falls under a category of effective field theory which is referred 
to as "rigorous potential models" in Ref. of which the most notable example is nonrelativistic QCD (NRQCD) 
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| p3| . NRQCD describes quarkonium states, which are heavy quark-heavy antiquark bound states. While the mass 
of each constituent is toq, the three-momentum of the relative motion is only of the order a s rriQ, and the kinetic 
energy is of an even lower order a 2 rriQ, where a s is the QCD coupling constant at scale itiq. Since a s is small, we 
have a separation of scales which allows us to expand the Hamiltonian in powers of a s . In particular, since the kinetic 
energy scale is much smaller than the three-momentum scale, it is natural to impose a stronger cutoff on energies than 
on three-momenta in the effective field theory. The resultant effective theory is local in time but not in space, i.e., 
a potential. Potential models constructed with such a philosophy are rigorous in the sense that they are related to 
the original field theory through Wilsonian renormalization. Our treatment of heavy baryons is similar to NRQCD, 
with A 1 / 2 playing the role of a 2 . The similarity is even more apparent when one realizes that our formalism, just like 
NRQCD, can be viewed as a nonrelativistic expansion. In NRQCD v 2 ~ a 2 , while in our case v 2 ~ A 1 / 2 . 

We will end on a cautionary note with a comparison with another "rigorous potential model" , namely the nucleon- 
nucleon effective field theory. The deuteron is a nonrelativistic bound state, with the binding energy ~ 2 MeV order of 
magnitude smaller than the nucleon mass ~ 1 GeV. However, the large N c counting rules would have suggested a very 
different picture. Since baryon-baryon interaction is of order N c , both the binding energy Vq and the spring constants 
k of a deut eron are of order N c . With reduced mass fi ~ N c , large N c scaling rules suggest that the excitation energy 
lu = k/ fi ~ AT?, which is much smaller than the binding energy Vq. This implies the existence of many bound states 
beneath the dissociation threshold (the number of bound states should be of order Vq/oj ~ N c ), and the low- lying 
bound states should be deeply bound. This is in blatant disagreement with the deuteron in the real world, which 
is barely bound with a tiny binding energy in comparison to Aqcd : Vq ~ 2 MeV in the triplet channel, and the 
singlet channel is not even bound. It turns out that there are many different physical contributions to the binding 
energy. Each of these contributions may be of order N c , but it happens that they almost cancel completely, and 
the numerical value for Vq turns out to be accidentally small. This illustrates a fundamental limitation of counting 
schemes: a physical quantity may carry a numerical value very different from what the formal power counting suggests 
due to accidental cancelation or appearances of unnaturally large or small coefficients. One should be aware of the 
possibilities of such accidents and carefully check if the physical picture suggested by the counting rules resembles the 
real world. 

For our analysis of the heavy baryon, the binding energy Vq is formally of order unity, while the excitation energy 
u> is of order A 1 / 2 . As a result, one expects the number of bound states to be of order X" 1 ^ 2 ~ N^ 2 — \f% in the 
real world. Experimentally Vo and u> have been determined to be about 625 MeV and 330 MeV, respectively. These 
numbers are at least compatible with the picture suggested in this paper. With an expansion parameter as large as 
A 1 / 2 - l/y/3 ~ 0.6, however, the expansion series may converge slowly and one needs to include corrections due to 
higher order terms before one can make any quantitative predictions. As a result, it is imperative to make a careful 
study of the higher order corrections. 

In summary, we have introduced an effective theory to study excited heavy baryons which makes the existence of 
the emergent symmetry manifest. What are the phenomenological implications of this emergent symmetry? What 
does this symmetry tell us about the strong decays of excited baryons and the weak decay form factors? And most 
importantly, are these symmetry predictions safe against corrections due to higher order terms in the A counting? All 
of these issues will be discussed in an upcoming paper |15[ ] . 
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